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Abstract 

We calculate the quantum states of regular polygons made of ID quantum wires treating each polygon vertex as a scatterer. The 
vertex scattering matrix is analytically obtained from the model of a circular bend of a given angle of a 2D nanowire. In the single 
mode limit the spectrum is classified in doublets of vanishing circulation, twofold split by the small vertex reflection, and singlets 
with circulation degeneracy. Simple analytic expressions of the energy eigenvalues are given. It is shown how each polygon is 
characterized by a specific spectrum. 
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1. Introduction 

Nanowires are a long-lasting topic of interest in nanoscience for two main reasons. One is, undoubtedly, their use 
as electron waveguides in devices and technological applications. The other is the possibility of using nanowires to 
artificially control quantum properties at the nanoscale, thus improving our fundamental understanding and prediction 
capabilities. In nanowires made with 2D electron gases of semiconductors like GaAs quantum behavior manifests at 
a mesoscomc scale, beyond the atomistic description, and it can be described with effective mass models and smooth 
potentials lHIS]- 

The quantum states on nanowire bends attracted much interest some years ago SlliiSSilil. It was 
shown that bound states form on the bend, and the transmission and reflection properties as a function of the energy 
and width of the wire were thoroughly investigated. In the single-mode limit of a 2D wire a bend can be described 
with an effective ID model containing an energy-dependent potential |@]. In the case of a circular bend, the situation 
becomes simpler and analytical approximations using a square well whose depth and length are fixed by the radius 
and angle of the bend were suggested ^ . 

Motivated by the above mentioned interest on nanowire bends we study in this work polygonal structures made 
of ID quantum wires (Fig. [T]). We focus on the single-mode limit and describe each vertex as a scatterer. The 
vertex scattering matrix is taken from the model of circular bends in 2D nanowires. We And that the spectrum fully 
characterizes the polygon structure. Two kinds of states are present: a) doublets, with a small energy splitting due 
to reflection and with a vanishing circulation along the perimeter of the polygon; and b) singlets, with an underlying 
degeneracy by which the circulation can take arbitrary values within a given range. The energy splitting of the doublets 
depends on the reflection probability of the vertex, while the number of singlets between consecutive doublets (in 
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Figure 1. a) Sketch of the regular polygons made of ID nanowires. The vertex labeling and the definition of a longitudinal coordinate x is indicated 
only for the triangle, b) Scattering amplitudes as defined for a given vertex i. 


energy) can be used to infer the number of vertices. The energy scale, setting the separation between states, is fixed 
by the length of the nanowire forming the sides of the polygon. 

We finally remark that, although our approach is on physical modeling of nanostructures, there is a connection 
with the more mathematically oriented study of quantum graphs 01 . Indeed, from this perspective our results can be 
viewed as a particular application to the case of regular polygons of the more general question of reconstructing the 
graph topology from the knowledge of the spectrum. 


2. Theoretical model and method 


We consider a regular polygon with vertices, made of ID quantum wires of length i (Fig. [T]). The distance 
along the perimeter is measured by variable x, with origin arbitrarily taken on the first vertex. The wave function 
between vertices i and / + 1 is a superposition of left- and right-ward propagating plane waves 




^ip(x-Xi+i) ^-ipix-Xi+i) 


( 1 ) 


where the first and second equalities take the reference point on vertices i and / -f 1, respectively. The and b)[ 


"//r ^llr 

coefficients are the characteristic input and output scattering amplitudes defined in Fig.[Tp. We are assuming a single 
propagating mode of wavenumber p (see [Appendix A| ) and the set of vertex positions is {x/, / = 1,..., Ay}. 

The i-th vertex scattering equation relates output and input amplitudes as 
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where t and r are the usual transmission and reflection complex coefficients. The scattering matrix in Eq. @ is 
summarizing the physical effect of the vertex by means of two complex quantities, t and r. As they are required 
inputs in our model, we will take these values from the known scattering matrices of circular bends in 2D nanowires. 
Actually, as discussed in [Appendix A[ in the single-mode limit there is an analytical description of the scattering by 
a circular bend in terms of an effective square well depending on the radius and angle of the bend isj] . 

There is a relation between input and output amplitudes of successive vertices. 
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The phase e in Eqs. © and Q appears due to the assumption of two different reference points when considering 
the scattering processes from two consecutive vertices. 
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Figure 2. Energy dependence of the T measure in arbitrary units (black) and the circulation C in units of hplm (gray-red). The unit of energy is 
Imd^. The wave number p and the wire parameter d are defined in Sec.[2| As indicated, each panel corresponds to a polygon with a different 
number of vertices Ny but the same arm length £ = 30d. The figure shows the presence of doublets with vanishing circulation and singlets with 
maximal positive circulation Co- Notice that at the energies of the singlets the circulation can take values from -Co to -i-Co. Other parameters: 
R= l.3d. 


We define the circulation 

c = (5) 

independent on the choice of vertex i due to flux conservation. The circulation is also the same if one chooses the right 
coefficients and instead of the left ones in Eq. ©. Physically, the circulation is measuring the probability 
current flowing along the polygon in a particular state, characterized by a set of a and h coefficients and a wave number 
p (corresponding to an energy E). 

In [Appendix B| we define an energy-dependent measure T{E), such that it vanishes for the eigenenergies of the 
polygon. In practice, we scan numerically the values of 7^{E) in order to determine the eigenenergies within a given 
energy interval. 


3. Results 


Figure [2] shows the energy dependence of T and C for polygonal structures from 3 to 6 vertices for a represen¬ 
tative set of parameters. As mentioned above, the energies for which vanishes are the physical eigenenergies of 
the polygon. We considered a representative energy interval above the threshold si = b?7i^I(2md^). The energy 
Si corresponds to the first transverse mode of the 2D nanowire of width d from which our ID model derives (see 


Appendix Aj ). 


From Fig. [2] we notice that the polygon spectra are characterized by a sequence of doublets, with small energy 
splittings and with vanishing circulation C. In addition to the doublets, there are a varying number of single-energy 
modes (singlets) with finite C lying in between doublets. For Ny = 3, for instance, doublets around 5.02 and 5.13 
with two intermediate singlets are seen in the upper left panel. Remarkably, doublets at the same energies are present 
in polygons with odd and even number of vertices. Polygons with an even number of vertices, however, possess 
additional doublets at intermediate energies. 

Singlets in Fig. [2] are characterized by having C ~ Nyhplm = Co- Notice, however, that different values of 
the circulation are possible for each singlet energy. Indeed, with the algorithm of Appendix Bj we find that a positive 
circulation Co is obtained if requiring = 1 (for an arbitrary i) and negative circulation -Co if requiring = 1. This 
indicates that for each singlet the circulation can actually take any value within the range -Co < C < Co by adequately 
superposing the solutions for positive and negative circulations. We have explicitly checked this possibility from the 
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Table 1. Different types of modes in vanishing-circulation doublets. The pairs forming doublets are (I,II) and (III,IV). 
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numerical solutions. This behavior manifests the connection between the symmetry breaking and the C-degeneracy 
of the singlets. 

A qualitative difference is seen when comparing left and right panels of Fig. [2l Polygons with an odd number 
of vertices have - I singlets in between vanishing-circulation doublets, while polygons with an even number of 
vertices only have 1. This explicitly shows that by knowing the spectrum one may characterize the polygonal 

structure. For instance, a spectrum consisting of a sequence of mode doublets with three intermediate singlets would 
correspond to an octagon. 

For the spectra with an even number of intermediate singlets a confusion might arise between the polygons with 
Ny and 2Ny vertices. For instance, both the triangle and the hexagon in Fig. [2| have two intermediate singlets. Still, 
one may differentiate the two situations with the discussion on doublets of next subsection. In essence, for odd-TVy 
polygons all doublets are similar in that each polygon side contains an even number of half wavelengths. On the 
contrary, for even-TVy polygons consecutive doublets alternate from an even to an odd number of half wavelengths. As 
shown below, this difference can be seen counting the number of density maxima on a side of the polygon. 

3.1. Doublets 

The states with vanishing circulation are characterized by having \a\ = \h\ on each side of the polygon. Besides 
vanishing circulation, doublets have a density that remains invariant by a translation from one polygon side to the next 
one (Figs. [3] and ID). This happens because the scattering amplitudes of successive vertices fulfill either or 

At the same time, for a given vertex the coefficients also fulfill ai = ±ar. Therefore, these vanishing- 
circulation states can be classified in four types as summarized in Tab. [T] 

Table [T] also gives the secular equation determining the wave number p and the energy of the eigenmode for each 
type of mode. Notice that solutions of types I and III are associated to states with maximal density on each vertex, 
while types II and IV have minimal density on the vertex, cf. Figs. [3] and (H The split pair forming a doublet are 
therefore composed by solutions of types (I,II) and (III,IV). It is also worth stressing that solutions of type III and IV 
are not allowed in polygons with an odd number of sides. The reason is most easily understood assuming r 0 and 
t ^ 1. In this case, the secular equation for types III and IV simplifies to = -1; that is ^ = {2n + l)/I/2 and an 
odd number of half wavelengths should fit in a polygon side. For odd values of Ny this is not compatible with the 
additional condition that the full perimeter should contain and integer number n' of wavelengths, Nyi = n'A. 

Simple analytical approximations for the doublet splittings can be obtained from the secular equations in Tab. [T] 
assuming t ^ I and r ^ This leads to the conditions (^ = 1,2, 3 ...) 

pi = 2nn + ro sin 0 , (I, II) , 

pi = {2n+\)7i + ro sin 0 , (III, IV) . 

For both kinds of pairs the splitting in wave numbers is Ap = 2ro sin 0/^, with an associated energy splitting 
A = /m)pAp. Equation (O clearly shows that in (I,II) doublets, the solution of type I decreases its energy due 
to reflection while that of type II increases, and analogously for (III,IV) doublets. We recall that type I and III solu¬ 
tions are those with a maximal density on each vertex. Finally, we notice that as the energy increases the reflection 
coefficient is reduced, and so the doublet splitting is also reduced, as hinted also in Fig.[2l 

Summarizing, the energies of the zero-circulation doublets are approximately 
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Physically, our result of Eq. © says that the arm length £ determines the separation between successive doublets, 
while the reflection properties of the vertex determine the splitting A of the pair forming a doublet. 


3.2. Singlets 

The most remarkable feature of singlets is that the density is not equivalent on different sides of the polygon. As 
mentioned above, for each singlet different solutions with circulation ranging from -Co to +Co, where Co ~ Nyhplm, 
are possible. The intermediate panels of Figs. [3] (energies of 5.049 and 5.086) and in Fig|4](5.042 and 5.097) show the 
densities corresponding to singlets for the triangle and square. Continuum and dashed lines are for states with maximal 
and vanishing circulation, respectively. We notice that for states with non vanishing C, the density for positive and 
negative circulations coincide. For nonvanishing circulation the density does not vanish at any point, but oscillates 
around a finite mean value. In all singlets, different behaviors are seen on the sides of the polygon. 

The singlet energies and their underlying C-degeneracy can be easily understood in the limit of vanishing reflec¬ 
tion. Indeed, when r 0 and 7 1 the solutions to Eqs. (O and Q decouple for positive and negative circulations. 

The input coefficients for vertex k read in this case 

a® = (C < 0), 

af = , (C > 0), ^ 


where a and [3 are arbitrary independent constants that can be fixed by normalization. A maximal circulation state is 
obtained when either a ox{3 vanish. Arbitrary superpositions can be formed for a 9 ^ 0 and (3 including a vanishing 
circulation state for a = JS. The cyclic condition in Eq. = 1 , leads to the energies 




(n= 1,2,3,...). 


(9) 


These energies coincide with a doublet condition lEq. d?]) for A = 0] each time n is an integer number of times Ny. It 
can be easily verified that Eq. (|9]) implies that for odd-Ay polygons there are Ny - I singlets between two successive 
doublets and that for even-Ay polygons there are Ny/2- 1, in agreement with the results of Fig. [2] 


4. Conclusions 

A model of closed polygons made of ID quantum wires has been presented where each vertex is described as a 
scatterer. The vertex transmission and reflection coefficients have been described with the model of circular bends 
in 2D wires. The polygon spectra are characterized by a sequence of doublets, with a small energy splitting, and 
a typical number of singlets lying in between doublets. Odd-Ay polygons have Ay - 1 singlets while even Ay have 
Ay/2 - 1 singlets in between two doublets. The doublet splittings are caused by the reflection on vertices. Doublets 
have a vanishing circulation and singlets can have circulations ranging from a negative to positive characteristic 
values. Approximate analytical expressions for the singlet and doublet energies have been given and the corresponding 
densities discussed. Our results explicitly show how the polygon characteristics can be inferred from the spectrum. 
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Appendix A. Scattering by a 2D circular bend 

We consider a 2D nanowire with a circular bend of angle 26 and radius R, as sketched in Fig. IA.5I Notice that 
6, defined as half the bend angle, is related in our model to the number of polygon vertices by ^ = n/Ny. The lateral 
extension of the nanowire is J, setting the energy of the first transverse mode to /(2md^). Following Ref. 

iH], it is possible to derive an analytical expression for the reflection and transmission coefficients of the circular bend 
in the single-mode limit. 
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Figure 3. Densities for the triangle modes of Fig.|3 All distances are measured in units of d (defined in [Appendix A) . Upper and lower panels 
are for vanishing-circulation doublets of type (I,II). Intermediate panels are singlets with maximal (solid) and vanishing (dashed) circulation. The 
energy of the mode (as in Fig. [3 is indicated in each panel. The position of vertices is given by the vertical lines. Rest of the parameters as in Fig. 

m 


Defining a = R6, where R = 5 /R(R + d) is an average radius, the 2D scattering problem in the single mode limit 
is approximated by a ID quantum well of width 2a and depth Vq = K^mR^). The analytical expressions of the t 

and r complex coefficients are given in terms of two wave numbers (E > si) 


Specifically, they read 


12m ^ ~ 

__ 

{p + ^)2e2'(P-?)a -(p- q)2e2iip+q)a ’ 

2i(p - q){p + q) sin(2^/a) 

(p + q)2(;2i{p-q)a - (p - q)2g2i(p+q)a 


(A.1) 

(A.2) 


(A.3) 

(A.4) 


Notice that, as expected, for large energies itis p ^ q and Eqs. (IA.3I) and (IA.4b then yield t ^ 1 and r 0. We also 
stress that, due to symmetry, the transmission and refiection coefficients for left and right incidences are identical, 
t' = t, r' = r, leading to the scattering matrix of Eq. 


Appendix B. Algorithm 

Equations Q, (O and 0 are set up as a homogeneous linear system 


0 , 
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Figure A.5. Circular bend in a 2D nanowire. 


where we define the vector 




(B.2) 


with an analogous definition for b. Nontrivial solutions of Eq. (IB.II) occur only for energies such that the determinant 
of M vanishes, det(M) = 0. 

In practice, we arbitrarily pick one of the components of a, say a^\ and fix it to 1, transforming Eq. (IB. II) into the 
inhomogeneous problem 


M' 


a 

b 



(B.3) 


where 

= (0...010...0...0), (B.4) 

with the value 1 on the position of the chosen component a^\ Matrices M and M' are identical except for the row of 
the chosen component, which is diagonal for M'. Equation (IB.31) can be easily solved by standard numerical routines. 
The solution is then used to evaluate the following norm 


‘F = norm 




(B.5) 


Obviously, when ‘F = 0 the solution obtained from Eq. (IB.31) is actually a solution of the original problem (IB.lb . 
The validity of the method is guaranteed by the linearity, for when a solution of (IB.lb exists it can always be scaled 
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such that a chosen component is equal to 1, provided only that it does not vanish. A simple energy scan of T will 
signal the position of the eigenenergies as the nodes of the function. We have checked that, as expected, the nodes are 
independent of the chosen component and that they numerically coincide with the nodes of det(M) = 0. 
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